Limiting SUSY-QCD spectrum and its application for decays of superheavy particles 
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The supersymmetric generalization of the limiting and Gaussian QCD spectra is obtained. These 
spectra are valid for i < 1, when the main contribution to the parton cascade is given by gluons 
and gluinos. The derived spectra are applied to decaying superheavy particles with masses up to 
the GUT scale. These particles can be relics from the Big Bang or produced by topological defects 
and could give rise to the observed ultrahigh energy cosmic rays. General formulae for the fluxes of 
protons, photons and neutrinos due to decays of superheavy particles are obtained. 
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I. INTRODUCTION 



The spectra of hadrons produced in deep-inelastic scattering and e + e annihilation are formed due to QCD cas- 
cading of the partons. In the Leading Logarithmic Approximation (LLA) which takes into account ln(Q 2 ) terms this 
CN| \ cascade is described by the Gribov-Lipatov-Altarelli-Parisi-Dokshitzer (GLAPD) equation 0. The Modified Leading 
• Logarithmic Approximation (MLLA) takes into account additionally ln(x) terms, where x = k\\/k™ ax and fc|| is the 
longitudinal momentum of the produced hadron. Color coherence effects arc described in MLLA. Two approximate 
analytic solutions to the MLLA evolution equations have been obtained. These are the limiting spectrum p] and the 
Gaussian spectrum ||,|| , in which we include the distorted Gaussian spectrum . The limiting spectrum is the most 
accurate one among them. In fact, it describes well the experimental data at large x, too, and this is natural, though 
I/") ' accidental (for an explanation see Ref. ||). The limiting spectrum has a free normalization constant, Ku m , which 
cannot be calculated theoretically and is found from comparison with experimental data. This constant has to be 
considered as a basic parameter of the theory, and it can be used at all energies, where the physical assumptions, under 
which the limiting spectrum is derived, are valid. For detailed calculations of hadron spectra in e + e~ annihilation and 
comparison with experimental data see Ref. S. Up to energies of existing e + e~ colliders, y/s < 183 GeV, the limiting 
I spectrum and the distorted Gaussian spectrum describe well the available data. At large energies y/s > 1 TeV the 
production of supersymmetric particles might substantially change the QCD spectra. Apart from future experiments 
at LHC, supersymmetry (SUSY) might strongly reveal itself in the decays of superheavy particles. They can appear 
as relics of the Big Bang, or be produced by topological defects (TD), and can be the sources of the observed ultrahigh 
' energy cosmic rays (UHECR) at E > 1 • 10 10 GeV. The range of masses, mx, of interest for UHECR goes from the 
> GUT scale (m x ~ 10 16 GeV) or less down to m x ~ 10 12 - 10 14 GeV. 

In this paper we obtain the generalization of the limiting and Gaussian QCD spectra for the SUSY case. Although 
the influence of SUSY on, e.g., the evolution of parton distributions or the running coupling constant was considered 



' in many works in the 80's H, this is, to the best of our knowledge, the first time that fragmentation spectra of hadrons 
are examined for large yfs up to the GUT scale in SUSY-QCD. As application we use these spectra for calculations 
of the fluxes of UHECR produced in the decays of superheavy particles. 



II. LIMITING SPECTRUM IN SUSY-QCD 

The GLAPD equation jl[] describes in LLA the evolution of the parton distributions D^(x,£,) with £. Here is 
the distribution of partons B inside the parton A dressed by QCD interactions with coupling constant a(k\), where 
k± is the transverse momentum and x — fc||/fc" lax is the longitudinal momentum fraction of the parton B. The 

variable £ characterizes the maximum value of k\ available in the considered process (k\ < Q 2 ), 




with A ~ 0.25 GeV as phenomenological parameter. 

In LLA, when terms with a s (Q 2 ) ln(Q 2 ) are kept and terms proportional to a(Q 2 ) are neglected, the GLAPD 
equation can be written as (|, Eq. (1.79)] 
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where & B (z) is the splitting function characterizing the decay A — > B + C . 

The supersymmetrization of Eq. (0) is simple: each parton A should be substituted by the supermultiplet which 
contains A and its superpartner A. We shall generalize here the limiting spectrum of QCD given e.g. in [0] to SUSY- 
QCD. This spectrum is valid for small x, where gluons strongly dominate the other partons. Therefore, we restrict 
ourselves to calculations taking into account only two partons, namely gluons g and gluinos A, in the tree diagrams. 
However, in the loop diagrams which govern the running of a s {k 2 _) we take into account also quarks and squarks. 

Multiplying Eq. (H) by x^ +1 and integrating it over x, we obtain an equation for the moments D B (j,£), 
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where 
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and the indices A, B, C run through g and A. 
The splitting functions are given e.g. in Q as 
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Equation (^) can be rewritten in matrix form choosing as basis {g, A), 
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where 
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The procedure of solving Eq. (||) is identical to that given in Ref. ||. After diagonalization of H(j) with the help of 
(D + ,D~), the eigenvalues of H are 



" = j \ v a + u * ± 
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In the limit uj = j — 1 — ► 0, the leading term v + is given by 
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with a — ^-N c = 11. (N c — 3 is the number of colours.) 

Up to now we have considered the LLA approximation. This approximation is not correct for x <C 1, when 
colour coherence effects become important. A better description of distributions at small x is given by the MLLA 
approximation ||. It takes into account both ln(Q 2 ) and ln(x) terms as well as angular ordering. The MLLA evolution 
equation results in Eq. (j^f for the moments, with £ replaced by £mlla = a s(Y)/(4:ir) ln(Y 2 ) with Y = Iyl(E9 /Qo), 
where E is the energy and 6 the opening angle of the jet (see || Eq. (7.4)]. For its analytic solution at small x, 
the limiting spectrum, the eigenvalues (13) derived above in LLA are still valid. Properly normalized, the limiting 
spectrum Dy lnl (l,Y) = xDu m (x,Y) gives a~ 1 da/dl in the case of e + e~ annihilation and the decay spectrum of X 
particles, where I = ln(l/a;). This spectrum is given by JlO[ 

DUl,Y) = K^T(B) r d -?- Ba {^— V y) B lB(y)- (15) 
J-n/2 T \&N C a YJ 

Here, Y now is Y — ln[ v / i/(2A)], a = ao + it, uq = arctanh(2£ — 1), £ = 1 — l/Y, and y/s is the cm. energy of an 
e + e~ pair or the mass mx of the superheavy decaying particle. The parameters depending on the structure of the 
theory are a = ^N c , the constant b of evolution of a s (k 2 J in one-loop approximation, k 2 ± da s {k 2 _) / dk\ = — &a 2 /(47r), 
and Cf = (N 2 — 1)/(2N C ) = 4/3. When the masses of the superheavy coloured Higgses are larger than Q 2 = k\ max , 
b = &susy = 9 — rif, where nt — 6 is the number of quark flavours. Finally, Ib is the modified Bessel function of 
order B = a/b, and argument 

y(r)=( ^£ [cosh a + (1 - 2Q sinh a] y\ . (16) 

\ b smha / 

A convenient way for the numerical evaluation of Ib{v) is the use of its series expansions, given for example by (8.445) 
and (8.451.5) in Ref. O). 

In Fig. 1 the SUSY-QCD limiting spectrum is shown in comparison with the QCD limiting spectrum, for the case 
of three colours, N c = 3, and six quarks flavours, nt = 6. For the sake of comparison, we normalize both spectra by 
the condition 

l 

dx xD nm (x, Y) = 2 . (17) 



The maximum of the spectra is at l m = Y~(0.5 + ^JTJy - c/Y) with c w 0.39 in QCD and c w 0.84 in SUSY-QCD. 
Therefore, the SUSY spectra are shifted to the right, and since they are also narrower than the QCD spectra (see 
Eq. (|l9| ) below), the SUSY maxima are dramatically higher (by a factor of 30) than the QCD ones. We remind the 
reader that the value of the maximum is given by the multiplicity. 

We compare these spectra also with the Gaussian approximation obtained for the MLLA solution HQ. This 
approximation can be easily generalized to the SUSY case |l^| and has for x -C 1 as function of x = 2E/^/s the form, 

K G ( ln 2 x/x„ 
— 6XP { 2^ 



Z? G (x,Y) = ^exp(-^f^i) (18) 



Assuming one- loop SUSY evolution of a s (k 2 _) with b = &susy and A = Aqcd, one has 



and 



A 



1/2 



(20) 



Since &susy = 3 (for rif — 6) is less than 6qcd = 7 (for the same rif), the SUSY-QCD peak is narrower than the 
QCD one. 

The Gaussian spectrum given by Eq. (|l8|) has its maximum at x = x m . It gives a less precise description than the 
limiting spectrum. One might expect that at large Y , when the higher momenta, skewness s and kurtosis k, become 
small (s ~ y~ 3 / 2 anc l K ^ y -1 / 2 ) the agreement improves. We have found that both for the case of ordinary QCD 
and SUSY-QCD the shape of the spectra differs substantially in the interesting range 10~ 6 < x < 10~ 2 (Fig. 2). 
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III. APPLICATIONS 



We shall apply now the limiting spectrum in SUSY-QCD to the calculation of the spectrum of ultrahigh energy 
cosmic rays (UHECR) generated by the decay of superheavy particles with masses 10 12 — 10 15 GeV. 

Let us discuss first the problem of the normalization of spectrum. We remind the reader that the limiting spectrum 
has been derived for x <C 1, though at least for small s it describes well the experimental data for x up to x ~ 1 
(sec H for discussion). In Ref. fiofl , the normalization constant was fixed by a comparison with experimental data on 
e + e~ annihilation to K\\ m — 2.6. This value is fixed by Dn m (l,Y) at maximum (see Fig. 1), i.e. by multiplicity. 

Since the shape of the spectrum and the position of its peak change dramatically, if one goes from QCD and 
y/s ~ 100 GeV to SUSY-QCD and y/s ~ 10 12 — 10 15 GeV, we cannot use this value of Kn m . Instead we use as 
normalization condition Eq. (^) replacing the factor 2 by 2fi, where i runs through N (all nucleons) and ir^ and 
7T° (charged and neutral pions), while /, is the fraction of energy carried by the hadron i. Note that the main 
contribution to the integral in Eq. ([l7]) comes from large values x ~ 1, where the limiting spectrum might have large 
uncertainties. However, this is, in our opinion, the most physical way of normalization. The numerical values of fi 
are unknown at large s. One can assume that f T w 1 — /lsPj where /lsp is the energy fraction taken away by the 
lightest supersymmetric particle (LSP). According to the simplified Monte-Carlo simulation of Q], /lsp ~ 0.4. For 
the ratio /jv//tt we use ~ 0.05 inspired by Z° decay. 

Let us assume that the decay rate of X particles fix in the extragalactic space does not depend on distance and 
time. Then taking into account the energy losses of UHE protons and the absorption of UHE photons due to pair 
production (7 + 7 — ► e + + e~) on the radio and microwave background, the diffuse flux of UHE protons and antiprotons 
is 

W*0 = ^ ^ r dt g D N (x g ,Y)^§^l , (21) 
2vr m x Jo dE 

where E g (E, t g ) is the energy at generation time t g of a proton which has at present the energy E an d x Q = 2E g /mx- 
Denoting the proton energy looses on microwave radiation by dE/dt — b(E, z), dE g /dE is given by 15 



^§^ = (l + , ff )ex P 



dz_ h , .1/2 f db(E,0) 
H (L + Z) { dE 



E=E g (z) 



(22) 



where Hq is the Hubble constant and z the redshift. In the case that the energy looses on microwave radiation are 
much larger than the adiabatic ones, Eq. ( |2^ ) reduces to 

dE g (E,z g ) ^ b(E g ,0) 

dE ~ b(E,0) ' 1 ' 

The diffuse spectrum of UHE photons can be calculated as 

7 7 (£) = — A 7 (£)— / —D w o(x,Y), (24) 
Ti" rn x J2E/ m ^ x 

where \-y(E) is the absorption length of a photon. For the numerical evaluation of I 7 , we use A 7 from Jig] . 
The neutrino flux depends generally on the evolution of the sources, 



hx(t) =M*o) ( j J ■ (25) 

The number of neutrinos with energy E produced per decay of one X particle can be approximately calculated as 

12 P 1 dr 

N v {E)n / —D^(x,Y). (26) 



Combining Eqs. (E5I) and (|26|), we obtain for the diffuse neutrino flux 



Iv{E)= 3^0) f z - d z(1 + z)h f <t D%±{x , Y ), (27) 
irH m x Jo J±m+z) X 



where 1 + z max (E) w m x /(AE). 



4 



In Fig. 3, the spectra of UHE protons, photons, and neutrinos are shown together with experimental data |l7| for 
the model of cosmic necklaces - topological defects which consists of monopoles connected by strings |1| . We use the 
model of Ref . [|l3| , where 

"x = < 28 > 

with r 2 /i = 5 • 10 27 GeV 2 and rtix = 1 ■ 10 14 GeV. The proton flux is suppressed at the highest energies as compared 
with the calculations of Ref. jL3|, where the Gaussian SUSY-QCD spectrum was used. The model corresponds to 
p = 0, or effectively to the absence of evolution for the integral over z in Eq. (p7|). 

In conclusion, we have calculated the SUSY-QCD limiting spectrum of partons in a jet. This spectrum considerably 
differs from that of ordinary QCD: the maximum of the Gaussian peak is shifted towards smaller x and the peak is 
narrower and higher. The limiting spectrum has been applied for calculations of the UHE fluxes of protons, photons 
and neutrinos from decays of superheavy particles in the Universe. The fluxes of these particles are different from 
those calculated with the ordinary QCD parton spectra. 
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I = ln(l/x) 



Figure 1: Limiting spectrum D]i m (l,Y) as function of I = \n(l/x) for SUSY- 
QCD (solid lines) and QCD (dashed lines), both cases for mx = 10 12 GeV 
(bottom), m x = 10 13 GeV (middle) and m x = 10 14 GeV (top). The QCD 
spectrum is scaled up by a factor 30. 




x = 2E/m x 



Figure 2: Ratio Aim (a, Y)/D G (x, Y) as function of a; = 2E/m x for SUSY-QCD 
(solid lines) and ordinary QCD (dashed lines), both cases for rax = 10 12 GeV 
and m x = 10 14 GeV. 
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Figure 3: Predicted fluxes from cosmic necklaces with r 2 /x = 5-10 27 GeV: proton 
flux (I p ), upper (/ 7 ,i) and lower (/ 7 ,2) limit for the gamma ray flux and neutrino 
flux (I v ). The data points are based on the compilation of Ref. [13]. 



